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Recently it was found, in a broad class of models, that the dark energy density may change its sign 
during the evolution of the universe. This may lead to a global collapse of the universe within the 



time tc 



- 10 years. Our goal is to find what bounds on the future lifetime of the universe can 



be placed by the next generation of cosmological observations. As an example, we investigate the 
simplest model of dark energy with a linear potential V{4)) — Vo{l + acf>). This model can describe 
the present stage of acceleration of the universe if a is small enough. However, eventually the field 
rolls down, V{(j}) becomes negative, and the universe collapses. The existing observational data 
indicate that the universe described by this model will collapse not earlier than tc ^ 10 billion years 
from the present moment. We show that the data from SNAP and Planck satellites may extend the 
bound on the "doomsday" time to tc ^ 40 billion years at the 95% confidence level. 
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I. INTRODUCTION 

Recent discovery of the acceleration of the universe 
P, 0, is often interpreted as a proof that the universe 
is going to expand forever. This is indeed the case if 
the acceleration occurs due to the existence of a positive 
vacuum energy density A (cosmological constant), con- 
stituting approximately 70% of the energy density of the 
universe po ~ 10~^^°M^ ~ 10~2^ g/cm^ today [3|. How- 
ever, the only existing theoretical model based on string 
theory and describing an accelerating universe in a state 
with a positive cosmological constant predicts that this 
state is metastable, and therefore acceleration of the uni- 
verse cannot persist for an indefinitely long time 5]. 

The situation becomes even more complicated in the 
models where the present acceleration of the universe is 
related not to the cosmological constant, but to a slowly 
changing energy density of a scalar field, called dark en- 
ergy. One of the first models of dark energy (and by far 
the simplest one) was proposed back in 1986 in Ref. [H, 
where it was suggested to replace the cosmological con- 
stant by the energy density of a slowly changing scalar 
field (p with the linear effective potential 



F(0) = yo(i + "<? 



(1) 



Here we use the units Mp = (SttG)^^/^ 1. If the slope 
of the potential is sufficiently small, aVo ^ 10^^^", the 
field (j) practically does not change during the last 10^° 
years, its kinetic energy is very small, so at least until 
the present stage of the evolution of the universe its to- 
tal potential energy V{(l)) acts nearly like a cosmological 
constant. The anomalous flatness of the effective poten- 
tial in this scenario is a standard feature of most of the 
models of dark energy 0, |1] ■ 



The main reason to introduce this model in Q was to 
address the cosmological constant problem. The main 
idea can be explained as follows. Even though the en- 
ergy density of the field in this model practically does 
not change at the present time, it changed substantially 
during inflation. Since </> is a massless field, it experi- 
enced quantum jumps with the amplitude H/2t: dur- 
ing each time H~^. These jumps move the field in 
all possible directions. In the context of the eternal in- 
flation scenario this implies that the fleld becomes ran- 
domized by quantum fluctuations: The universe becomes 
divided into an infinitely large number of exponentially 
large parts containing all possible values of the field (p. 
In other words, the universe becomes divided into an in- 
finitely large number of 'universes' with all possible val- 
ues of the effective cosmological constant A{(j)) = F(0). 
This quantity may range from —M^ to +Mp in differ- 
ent parts of the universe, but we can live only in the 
'universes' with |A| < O(10)po ~ 10"^^ g/cm^. 

Indeed, if A < —10^^^ g/cm^^, the universe collapses 
within the time much smaller than the present age of 
the universe ~ 10^° years 0,0,^3- On the other hand, 
if A ^ 10~^* g/cm'^, the universe at present would ex- 
pand exponentially fast, energy density of matter would 
be exponentially small, and life as we know it would be 
impossible 0, Q ■ This means that we can live only in 
those parts of the universe where the cosmological con- 
stant does not differ too much from its presently observed 
value |A| ~ po. This approach proposed in constituted 
the basis for many subsequent attempts to solve the cos- 
mological constant problem u sing the anthro pic principle 
in inflationary cosmology [II [11 [11 [11 111 111 111 13 • 

However, the simplest dark energy model with the lin- 
ear potential JQ) has a disturbing consequence: even 
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though the field </> moves down very slowly, eventually 
the potential energy density becomes negative, and 
the universe collapses, just as the universe with a nega- 
tive cosmological constant. A detailed description of this 
process can be found in jl^j. Nevertheless, since the col- 
lapse will occur only in a distant future, the linear model 
([T|) with a sufficiently small a is quite satisfactory from 
the point of view of all existing observational data. 

The existence of the vacuum instability leading to a 
global collapse of the universe is a property of a large class 
of the models of dark energy, & Z 2L 2Z, 22.. 2M- 
However, it takes some time for the universe to switch 
from acceleration to collapse. Recent observational 
data allow us to rule out many of the models with 
steep potentials predicting rapid collapse of the uni- 
verse. Most of these models compatible with the exist- 
ing observational data predict the global collapse of the 
universe within the time exceeding 10-30 billion years 
El 113, m m m m m. This is the time compa- 
rable with the present age of the universe Iq ~ 13.7 
billion years. In this respect it becomes very inter- 
esting to check whether one could find any indication 
of the vacuum instability and the future "doomsday," 
or, vice versa, whether it is possible to increase the 
"life expectancy" of the universe, using the most ad- 
vanced data to be obtained, for example, by the Super- 
nova/Acceleration Probe (SNAP: 0|) distance-redshift 
measurements of supernovae,SNAP[SN], the Planck Sur- 
veyor cosmic microwave background satellite p^ . and 
weak gravitational lensing from the SNAP wide field sur- 
vey [ig, SNAP[WL], or the LSST survey 

This is the main goal of our paper. We will study this 
issue in the context of the simplest dark energy model 
. The reason to do it is that this model (up to the field 
redefinition) has only one free parameter a, so one can 
relatively easily study this model in all possible regimes. 
Also, this model is quite t ypic al: in most of the models 
discussed in 'T^, [2ll il ll [11 113 the potential looks 
linear with respect to (j) near the point where V{(j)) = 0, 
where the acceleration gives way to a collapse.^ 



II. DARK ENERGY WITH A LINEAR 
POTENTIAL 



We assume that the scalar field </>, with the linear 
potential (^3), represents the dark energy density of the 
universe. There is also the usual matter energy density 



Pm = f^, so that the equations of motion are given by 



^ A future curvature singularity may also appear in the models 
where the null energy condition is violated, p + p < , so that 
w < —1. We will not discuss such "phantom" models 1301 here 
because their physical interpretation is rather obscure and they 
are expected to lead to a very rapid development of instability 
at the quantum level [3l(. 



a 
a 



a o<p 



(2) 



(3) 



Here a{t) is the scale factor of the flat FRW metric, ds^ 
dt^ — a{t)^dx'^. The Hubble parameter is given by 



PM{t) + pD{t) _ pT{t) 



(4) 



We will solve these equations numerically, more details 
on this can be found in [2^. The dark energy density 
Pd and the pressure po are given hy po = 4>^ /'^ + V and 
Pd ~ <\>^ 1'^ — y . and the total energy density includes 
also the energy density of matter, pr = PM + Pd- 

A dimensionless dark energy (matter) density is given 
by a ratio of the dark energy (matter) to the total energy. 



Pt 



M 



Pm 
Pt 



(5) 



Observations suggest that now fijv/ ~ 0.28, flo ~ 0.72, 
and ^It — + ~ 1- The present time will be 
specified in our numerical solutions by the moment when 
flo = 0.72. We will study separately the effect of chang- 
ing the current value of between 0.7 and 0.73. An- 
other important characteristic of the dark energy is its 
pressure-to-energy ratio defining the dark energy equa- 
tion of state: 



WD 



PD 
PD 



02/2 - V 
02/2 + V 



(6) 



In what follows we will drop the index D in wd and use 
w for dark energy wu. 

Without any loss of generality one can assume that the 
initial value of the field (j) in the linear potential is zero, 
since any change of 4>q can be absorbed into a redefinition 
of Vq. The value of the constant part of the potential Vq 
for any choice of the slope is not independent: it is chosen 
in a way that today at z = the value of fljj equals 
0.72. Thus there is only one independent parameter in 
the linear potential model, a, or, equivalently, the slope 
of the potential, aVo. We will assume that a > 0, but 
all results for the lifetime of the universe depend only on 

We solved the equation of motion of the theory nu- 
merically and we present below a set of solutions of the 
linear potential model with various slopes. The scale 
factor a{t) and the equations of state w{z) are plotted in 
Figs, n and |21 respectively. These models are designed 
for the analysis of the current and future observational 
data. Complete information on each model is given in 
the TableQ] This includes the color of the curve, the val- 
ues of a and Vq, the slope aVo, the value of w at present 
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TABLE I: Parameters for the Models in Figs.[T]andH 



Parameter 


Cosmological 
Constant 


SNAP[SN] 
+ Planck 
+ SNAP[WL] 
( 9570 cl j 


SNAP[SN] 
+Planck 
( 95% cl ) 


SNAP[SN] 
( 95% cl ) 


Minimum 
Lifetime Model 


curve color 


red 


orange 


purple 


blue 


black 


a 





0.71 


0.76 


0.86 


1.13 


Vo 


0.72 po 


0.83 Po 


0.85 Po 


0.91 Po 


1.77 Po 


aVb 





0.72 X 10"'^" Af^ 


0.79 X 10~'^"Mp' 


0.96 X 10"'^" 


2.46 X 10''^" 




-1 


-0.89 


-0.87 


-0.82 


-0.0001 


tc 


oo 


39.5 Gyr 


35.5 Gyr 


28.7 Gyr 


11.3 Gyr 




t 



FIG. 1; Scale factor a{t) in five models, the present mo- 
ment is t = 0. The upper (red) curve corresponds to the 
cosmological constant model with the vanishing slope a; clas- 
sically it has an infinite future lifetime. The curves below 
(orange, purple, blue, and black) correspond to a steepening 
slope. The time remaining from today to the future collapse 
in these models is shown in the table. Time is given in units of 
Hg^ ^ 13.7/0.983 billions of years. In these units the current 
age of the universe t « 13.7 Gyr is given by 0.983. 




FIG. 2: Evolution of dark energy equation of state w{z) in 
five models; the present moment is at z = 0. The value w{0) 
is also given in the table. 



at z = 0, and, finally, the time tc of the collapse of the 
universe (from the present moment). 

Model 1 , the curve with ever expanding universe in Fig. 
n is a fiducial cosmological constant model with vanish- 
ing slope a; it has an infinite future lifetime. Models 2, 3, 

4 with increasing slope will later be associated with some 
limits on lifetime based on specific observations. Model 

5 has the largest slope aVo = 2.46 x 10~^^°Mp, for which 
the value of flo — 0.72 is barely reached. Any further 
increase of the slope will make the model ruled out by 
the data (assuming that ~ 0.72 at present) since 

in these models will never reach 0.72. The scale factor 
of the minimal model is plotted in Fig. 1, where one can 
see that the universe will collapse in a time of the order 
tc = 11 Gyr from now. 

tc 

35 
30 
25 
20 
15 
10 
5 

a 



0.2 



0.4 



. 6 



0.8 



FIG. 3: Dependence of the lifetime of the universe (starting 
from the present moment), in units of H^^) on a in the linear 
model. 

Fig. 121 shows the dependence of the time tc of the col- 
lapse of the universe (from the present moment) on a. 
As we see, the lifetime sharply increases for a <C 1. 

By looking at Figs. n][21one can easily conclude that 
the models with the lifetime smaller than 10 billion years 
are at odds with the existing observational data. First of 
all, in these models one would have flo < 0.72. Secondly, 
while current data cannot see time variation in w{z), the 

J da Oi5(a) w(a) 
J da (a) ' 



average value of w, defined as w 



would 



4 



be rather large, w > —0.6. However, for less extreme 
models, the average value of w is generically quite close 
to —1, so data on this quantity is not a reliable guide. 
One requires an experiment that is capable of seeing the 
time variation in w{z). This requires a new generation 
of surveys, which we discuss in the next section. 

Instead of trying to find out the best constraint on the 
lifetime of the universe on the basis of the present obser- 
vational data, we will try to understand to what extent 
the future data can allow us to predict the fate of the 
universe. Our results will be limited to the simplest lin- 
ear model, but they will be quite indicative of the general 
situation. 



III. CONSTRAINTS ON THE LIFETIME OF 
THE UNIVERSE FROM NEXT GENERATION 
OBSERVATIONS 

To constrain the dark energy models and their impact 
on the fate of the universe we use observational data from 
next generation cosmological probes. The centerpiece is 
the Supernova/ Acceleration Probe (SNAP) satellite, pro- 
posed as a dedicated dark energy mission designed to 
measure precise luminosity distances to some 2000 Type 
la supernovae covering the redshift range z = 0.1 — 1.7. 
This data is supplemented by another part of the primary 
mission, involving wide area measurements of weak grav- 
itational lensing. Valuable complementarity is provided 
by data from the Planck Surveyor cosmic microwave 
background satellite that will use the temperature fluctu- 
ation power spectrum to make precision determinations 
of several combinations of cosmological parameters. For 
our purposes the key quantity will be measurement of 
the angular diameter distance to the CMB last scatter- 
ing surface at z = 1089. 

Of course, we do not know the results of these future 
experiments. They may indicate that the universe al- 
ready began decelerating, which could be a precursor for 
the future collapse of the universe. Here we would like 
to study the most optimistic possibility. Let us assume 
that the new observational data will favor the simplest 
of all possible options: the universe is dominated by the 
positive cosmological constant with the equation of state 
w = —1. In terms of our simplest model of dark energy 
with a linear potential this would mean that the slope of 
the potential cannot be distinguished from zero with the 
accuracy of the combined set of experiments mentioned 
above. What kind of predictions for the future evolution 
of the universe we would be able to make? In particular, 
we wonder whether we will be able to say, as is often 
asserted, that the acceleration of the observable part of 
the universe will continue forever, and in about 150 bil- 
lion years our galaxy and its closest neighbors will remain 
the only inhabitants of the otherwise empty observable 
part of the universe. 

In order to study this question, we examine the types of 
cosmological information we will obtain. The luminosity 



and angular diameter distances are both related to the 
comoving distance by factors of 1 + z. The comoving 
distance is given in terms of quantities from Section ^ 
by 

The Friedmann equations of general relativity define 
the expansion history of the universe, the evolution of 
the scale factor a{t), in terms of the components of the 
energy density by 

{a/af = H^ = H^[nlj{l + zf + noiz)] (8) 

The dark energy density evolves with redshift z = 

a^^ — 1 as 

iloiz) = f7?,e3/c:['^^'/(i+^')l[i+"'(^')l (9) 

^ f70^(l + 2)3(l+"'o+«.a)g-3«,„^/(l+z) (^Q) 

The second line gives the result for a commonly used 
parametrization w{z) — wq + Wa (1 — a) giving a good 
approximation for slowly rolling fields [sj. Here Hq, ^^j, 
and are the values of H{z), Um{z), and ^Id{z) today 
at 2; = 0, respectively. 

Since the equation of state enters via two integrals, 
small deviations between the model behavior and the 
WQ — Wa fit will be unimportant. This allows us to plot re- 
sults in the wq — Wa space, easing comparison with other 
dark energy models. Note that a measure of the time 
variation of the dark energy equation of state is often 
given by w' = dw/d\n{l + z)\z=i = Wa/'i- 

Specifically, we take the SNAP baseline mission pre- 
sented in [33I, including statistical and systematic errors 
amounting to 1% in distance at the depth of the survey, 
z = 1.7. We marginalize over the absolute magnitude pa- 
rameter M. (including the Hubble constant i/o) and take 
a fiducial cosmological constant model with 17 = 0.72, 
see Appendix for details. For those cases where we only 
consider supernova data, we also impose a gaussian prior 
on VId of 0.03; when we include CMB or weak lensing 
(WL) data, the natural determination of Sl^ in comple- 
mentarity with the supernovae is better than this. When 
we incorporate WL data we use only information from 
the linear part of the mass power spectrum, in the man- 
ner of 34] . Note that WL is roughly equivalent to the 
CMB in complementary power with the supernova data. 
But if the fiducial model has time varying equation of 
state 'w{z) then there is a further gain in precision with 
all three data sets. 

Constraints on the dark energy model from the data 
are analysed within the Fisher matrix method j85j | which 
is also explained in the Appendix. The equation of state 
function w{z) for each model is represented by 2 param- 
eters, Wq and Wa- The fit suggested in js^ . 

W{z) ^ WQ+Wa{l- a) ^ Wq + Wa— , (H) 
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TABLE II: Parameters for the Ellipses and Models in Fig. El 





Cosmological 
Constant 


SNAP[SN] + 
Planck + 
SNAP[WL] 
( 68% cl ) 


SNAP[SN] + 
Planck + 
SNAP[WL] 
( 95% cl ) 


SNAP[SN] + 
Planck 
( 68% cl ) 


SNAP[SN] + 
Planck 
( 95% cl ) 


SNAP[SN] + 
( 68% cl ) 


SNAP[SN] + 
( 95% cl ) 


curve 
color 


red 


orange 
dashed 


orange 


purple 
dashed 


purple 


blue 
dashed 


blue 


a 





0.576 


0.71 


0.63 


0.76 


0.72 


0.86 




0.72 pn 


0.79 po 


0.83 po 


0.80 po 


0.85 po 


0.84 po 


0.91 Po 


10(0) 


-1 


-0.94 


-0.89 


-0.92 


-0.87 


-0.89 


-0.82 




CO 


55.3 Gyr 


39.5 Gyr 


47.9 Gyr 


35.5 Gyr 


38.6 Gyr 


28.7 Gyr 



works well for our models, especially for the models with 
larger lifetime. Plots in the wq — Wa plane marginalize 
over the value oi^o- 

Each case that we studied is shown in a wq — Wa plane 
as a (black) triangle in Fig.^] Each case (apart from the 
cosmological constant) is chosen so that the point in the 
Wo — Wa plane is at the boundary of one of the six ellipses, 
corresponding to future data. Complete information on 
each case is given in Tables I and ^ The color code 
for the confidence ellipses, orange, purple and blue, is 
the same as in previous figures, for SNAP-|-Planck-|-WL, 
SNAP+Planck, and SNAP, respectively with 95% confi- 
dence. Dashed lines in the same colors are used for 68% 
confidence data. We give the values of a, Vb and wq, and 
the lifetime before the collapse for all these cases. 

The triangles at the boundary of each ellipse give us 
information about the constraint on the lifetime before 
collapse in each case. Any model whose point in the 
Wo — Wa plane lies outside the corresponding ellipse will 
be ruled out if the actual observation will favor the cos- 
mological constant as the most likely point in the pa- 
rameter space. This will rule out models with lifetimes 
before the collapse smaller than that of each model at the 
boundary. When drawing the ellipses around the (future, 
expected) measurement point, we made the common as- 
sumption that the probability distribution has the same 
shape around the measurement value as around the true 
value. This assumption is justified, as, indeed, we con- 
firmed that ellipses drawn around the triangles are very 
similar in shape and size. 

The ellipses in wq — Wa parameter space for SNAP, for 
SNAP and Planck, as well as for SNAP and Planck and 
WL given in Fig. ^ are generic, model-independent, and 
can be used in connection with any model — not just the 
linear potential studied here — that has its w{z) reason- 
ably approximated by (|11() . The only information that 
enters are the properties of the measurements and the 
parametrization of the fit for w{z) used. Therefore our 
results can be easily used for investigation of other mod- 
els of dark energy. 

The lifetimes, on the other hand, that become associ- 
ated with each point in the wq — Wa parameter space de- 
pend on the individual model studied. The Fisher ellipses 
for SNAP, SNAP+Planck, and SNAP-|-Planck-hWL are 
shown in Fig. 0] Thus, the lifetimes that we attribute to 



the triangles there are particular for the linear model. 



Wa 

1 . 5 r 




-1 . 5 



FIG. 4: Confidence contours for different combinations of 
data sets are plotted assuming a fiducial cosmological con- 
stant (mo = —1, Wa = 0) model. The innermost, orange 
(dashed) ellipse represents SNAP[SN] + Planck + SNAP[WL] 
at 95% (68%) confidence level. The slightly wider purple el- 
lipses use only SNAP[SN] -I- Planck, and the rounder, blue 
ellipses use only SNAP[SN]. 

Thus our results can be formulated as follows. If the 
supernova distance-redshift observations from SNAP will 
yield the positive cosmological constant, i.e. wq = — 1, 
Wa = 0, as the most likely point in the WQ — Wa parameter 
space, our analysis allows us to rule out a future lifetime 
before collapse that is shorter than 28.7 Gyr at the 95% 
confidence level (solid blue ellipse), and 38.6 Gyr at the 
68% confidence level (dashed blue ellipse). 

With the future Planck mission added to the SNAP 
data, the corresponding lifetime constraints can be raised 
even further to 35.5 Gyr (sohd purple) and 47.9 Gyr 
(dashed purple) at the 95% and 68% confidence levels, 
respectively. 

Weak lensing similarly tightens the constraints on 
w{z). Both WL and Planck successfully complement 
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SNAP (but not each other). The three data sets together 
lead to 39.5 Gyr (sohd orange) and 55.3 Gyr (dashed or- 
ange) at the 95% and 68% confidence levels, respectively. 



IV. DIFFERENT fi?. 

While all our models include as a parameter to 
marginalize over, the fiducial, central value so far has 
been fi^j = 0.72 as suggested by a combination of CMB 
and large scale structure results p) and new supernova 
data 

It is interesting to estimate what the difference will 
be in lifetimes if one changes the fiducial to 0.7 or 
0.73. We have constructed the Fisher ellipses for SNAP 
and SNAP -f Planck for this case, see Fig. [S| We have 
also evaluated the relevant change in the lifetime bounds 
between cases fi^j = 0.7 and = 0.73. The bound is 
changed by few percent, therefore the changes in H.'jj do 
not lead to significant changes in expected bounds on the 
lifetimes. 




-1.5 



FIG. 5: Saturn plot of Fisher ellipses: Each of the Fisher 
ellipses shown in Fig. I^is now given for fiducial Q% taking 
3 values: 0.70,0.72,0.73. The innermost ellipse in each case 
corresponds to 0,% = 0.73, the outermost one corresponds to 
n°B = 0.70. 



V. POSSIBLE SIGNATURES OF THE FUTURE 
COLLAPSE 

Until now, we considered the possibility that the future 
observations will produce the simplest result, wq = — 1, 
Wa = 0, which would suggest that the dark energy is 
nothing but the cosmological constant Vq. In this case. 



because of the observational uncertainties, we would be 
unable to claim that the universe is going to accelerate 
forever, but we will be able to say, that in the context 
of the simplest dark energy model (1) it will not collapse 
earlier than in 40 billion years from now, at the 95% 
confidence level. 

1.5r 



1 




FIG. 6: Fisher ellipses with the centers corresponding to three 
different possible results of the future observations, including 
SNAP[SN] + Planck + SNAP[WL], see the text. The points 
outside the ellipses are ruled out at the 95% confidence level. 
The point wo = —1, Wa = 0, corresponding to the cosmologi- 
cal constant, lies outside of the green and blue ellipses and at 
the boundary of the red ellipse. 

But what if we get a different result? Suppose, for ex- 
ample, that the future observations (including SNAP[SN] 
+ Planck -I- SNAP[WL]) will tell us that wq = -0.83 and 
Wa = —0.26, which corresponds to our linear model with 
a — 0.86. This would imply that at the 95% confidence 
level the true values of wq and Wa lie inside the green 
ellipse with the center at wq = —0.83, Wa = —0.26 (see 
Fig. EJ. The blue ellipse corresponds to the possibil- 
ity that Wo = —0.88, Wa — —0.18, and the red one to 
Wo — —0.93, Wa — —0.1. 

Note that the point wq = —1, Wa = 0, corresponding 
to the cosmological constant, lies outside the green and 
blue ellipses and at the boundary of the red ellipse. 

This result has two different implications. First of all, 
if the future observations find any of the sets of wq and 
Wa discussed above, they will rule out the standard cos- 
mological constant model at the 95% confidence level (in 
the Fisher matrix approximation). 

Secondly, within the simplest model of dark energy, 
which is our linear model (1) P|, finding, e.g., that 
Wo = —0.93 and Wa = —0.1, would imply, at the 95% 
confidence level, that our universe will not exist forever, 
but is going to collapse. This would not be a definite 
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proof of the coming collapse, because there may be other, 
more complicated models of dark energy with similar val- 
ues of Wo and Wa, which do not lead to a global collapse. 
Still, this would be a serious warning sign, which would 
stimulate further investigation of the fate of the universe. 



VI. DISCUSSION 

Our investigation leads us to the following set of con- 
clusions. 

First of all, now it becomes even more apparent that 
it is very difficult to predict the future evolution of the 
universe. The standard textbook illustrations showing 
that open and flat universes slow down but expand for- 
ever, and a closed universe collapses, recently became 
replaced by the picture of a flat eternally accelerating 
universe. Now we are coming to a realization that the 
stage of acceleration may be transient, and even a flat 
universe may experience global collapse within a time 
comparable with its present age. 

Our investigation shows that even the best experiments 
to be carried out in the next decade probably will be un- 
able to give us a final answer concerning the destiny of 
the universe. Even if all experiments will unambiguously 
support the simplest possibility that dark energy is noth- 
ing but a positive cosmological constant with w{z) = —1, 
this will not really mean, as often claimed in the popular 
press, that in 150 billion years our galaxy and its imme- 
diate neighborhood will remain the sole island of matter 
surrounded by eternally expanding empty space. For all 
we know now, and for all we are going to learn in the 
next ten years, we will be unable to rule out the possi- 
bility that our part of the universe is going to collapse in 
the distant future. 

But one can look at it from a different perspective. 
We live at the very beginning of the era of precision cos- 
mology. It is amazing that within the next decade, by 
combining several different tools such as investigation of 
supernovae, CMB and weak lensing, we will be able to 
learn quite a lot about the possible outcome of the uni- 
verse evolution. For example, if the combination of the 
experiments discussed in our paper will show that the 
most probable parameters of the dark energy correspond 
to the simplest cosmological constant scenario, Wq = —1, 
Wa = 0, we will be able to say that in accordance with 
the simplest theories of dark energy, such as our linear 
model, there is no imminent danger of the global collapse 
at least for the next 40 billion years. 

On the other hand, if these observations will favor a 
different set of parameters, e.g. wq > —0.93, Wa — —0.1, 
this will rule out, at the 95% confidence level, the sim- 
plest cosmological constant scenario. This result would 
be of great importance for our understanding of the most 
fundamental issues of physics, such as the structure of the 
vacuum state. In addition, such a result would mean, at 
least in the context of the simplest model of dark energy 
(Q, that our universe is going to collapse. This would 



make further investigation of dark energy even more ur- 
gent and interesting. 
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Appendix: SNAP, Planck, and Fisher ellipses 

This appendix is designed to offer a simple, practical 
guide on the implementation of the Fisher matrix method 
for analysis of data constraints on cosmological parame- 
ters. For the mathematical basis of the method see [s^l 
and for a general application to cosmology see [s^ • Here 
we give a step by step introduction to allow those unfa- 
miliar with the method to use it immediately (also see 
the appendix of 33]). 

Formally, the Fisher matrix is defined by the expecta- 
tion value 



d"^ lnL(x,p) 
dpidpj 
d In L(x, p) d In L(x, p) 



(12) 



where L(x, p) = n^j^/(a;fc, p) is the combined probabil- 
ity distribution, and /(a;fc,p) is the probability distribu- 
tion of the individual measurement Xk that in general also 
depends on all the model parameters p = (pi, . . . ,pn)- P 
denotes the fiducial or (unknown) true parameter value 
of p. The second equality in H12(l follows from the matrix 
properties in the maximum likelihood approach, see [3^ 
and references therein for a derivation. In many cases 
the probability distribution L could be approximated by 
a Gaussian and its relation to observed quantities, such 
as supernova magnitudes, becomes very simple 37] . 

Practically, the Fisher matrix method provides a local 
approximation to the likelihood surface for model param- 
eters Pi given a set of observations x — {mk}- It involves 
the sensitivities, or first derivatives, dmk/dpi evaluated 
at the fiducial model {pi} [3^1. This means that it only 
gives good estimates for small fit uncertainties, i.e. where 
the data constraints are strong enough to limit consider- 
ation to models near the input. Biases as well as uncer- 
tainties can be treated within the Fisher method. 

The sensitivities combine into the symmetric Fisher 
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matrix: 



J- 7,1 



E 



1 drrik drrik 
a^{mk) dpt dp J 



(13) 



Here we have assumed that the covariance of the data 
points vanishes, so the covariance matrix of the errors 
reduces to diagonal entries l/a'^(mk)- This is a reason- 
able approximation under certain circumstances, such as 
when the data is grouped into redshift bins wider than 
the correlation length. Remember that the Fisher matrix 
is a quick and rough approach; if one wants a more rig- 
orous treatment then one might as well use a full Monte 
Carlo of the data set (including covariances) . 

One can see that the larger the derivatives (i.e. more 
sensitivity), the larger the matrix entries. Similarly, the 
more (independent) data points or the smaller the data 
errors, the larger the entries. This leads to the Fisher 
matrix also being known as the information matrix and 
this name gives a good guide to the interpretation. The 
larger a matrix entry, the more information the data has 
provided on those model parameters: hence it is a more 
sensitive probe and a better final parameter estimator. 

The parameter estimations and their covariances can 
be read off from the inverse of the Fisher matrix, known 
as the error or covariance matrix C — F~^. So for exam- 
ple Cii = crf^ gives the one sigma error (68% confidence 
level) on the parameter pi, averaged over all the other 
dimensions of the parameter space. Formally, the Fisher 
method only gives a lower limit on this error, as given by 
the Cramer-Rao inequality 



Ap, > 



1 



(14) 



known as marginalization, excises the appropriate rows 
and columns from the covariance matrix. Thus, as stated 
before, Cu gives the error on a single parameter, averag- 
ing over all others. 

If one has multiple data sets, then the information sim- 
ply adds (if the data sets are independent). That is, one 
simply adds the Fisher matrices, assuming they have a 
common parameter set (note, the m's in l|13|l can be com- 
pletely different quantities for the two Fisher matrices 
added, what is important is that the piS have exactly the 
same meaning) . So if one wants to constrain cosmological 
parameters using both supernova distance data and CMB 
power spectrum data, then one just adds the Fisher ma- 
trices of each experiment. As a shortcut, sometimes one 
directly places a parameter constraint on the phase space 
rather than incorporating the data from the beginning. 
This is known as a prior, and is also implemented by 
adding a Fisher matrix; if the prior is on a single param- 
eter, e.g. ill) is known to coj^, then the addition is of a 
matrix empty except for one entry Ff^^°'' — 1 /<Jq^ , where 
i = Strictly speaking, this should be done only if the 
prior knowledge truly has no covariance with any other 
parameters. If one is interested in the effects of different 
priors on the final parameter estimation, then one can 
use rules of matrix algebra to quickly determine how this 
extra bit of information affects the results (see js^l)- Bi- 
ases, where the data has been skewed from the true model 
behavior by systematic errors, i.e. ruk rnu + 5mk, can 
also be treated by matrix manipulation (see |39j | and the 
Appendix of |3^1. 



SNAP 



where Api = \/ {pf} — {Pi}^ is the variance of the 
model parameter p^, given the measurement errors. The 
Cramer-Rao lower bound theorem H14|l is proved readily 
from the general property 



Var(A:) > 



(Cov(X, Y))^ 
Var(y) 



(15) 



which holds for any two random variables X and Y, by 

taking X ^ P,{x,,...,Xk), Y ^ 91n(L(.,......^;p,....,p.)) ^ 

and using that, for these X and Y , we have (Y) — 0, 
Cov(A:, Y) = 1 and 



Var(r) 



d"^ ln(L(xi, . . .,xk]P1j ■ ■ ■ ,Pn)) 
dpi 



(16) 



One can deal with the full set of parameters in two 
basic ways: by fixing those not of immediate interest 
or by averaging over their probability distribution. The 
first is accomplished practically by excising the rows and 
columns corresponding to the fixed parameters from the 
Fisher matrix before inverting it to find the error ma- 
trix. This is a severe step, equivalent to assuming per- 
fect knowledge of that variable. The second approach, 



The SNAP data set will comprise high precision and 
accuracy measurements of the astronomical magnitudes 
of some 2000 Type la supernovae from redshift z = Q 
out to redshift z = 1.7. Magnitudes are logarithmic dis- 
tance variables and depend on the cosmological variables, 
say flo, wq, and Wa, and an unimportant (for cosmol- 
ogy) absolute supernova luminosity variable written as 
M . We will always marginalize over M , assuming a uni- 
form probability distribution. 

For the level of statistical analysis considered here, the 
supernova magnitude data is placed into 17 redshift bins, 
each bin having a width of (5z = 0.1. This should be wide 
enough that observational uncertainties are uncorrelated 
from bin to bin. The parameter analysis turns out to 
be not very sensitive to the exact number of supernovae 
within each bin, as lon g a s the entire redshift range is 
reasonably represented liol . The SNAP baseline distri- 
bution is tabulated in |33|: generically one also includes 
data from 300 expected supernovae of the currently run- 
ning Nearby Supernova Factory |4]| , as a single point in 
the lowest bin. 

The magnitude error a(mk) comprises two contribu- 
tions: a statistical component from observational and in- 
trinsic supernova magnitude dispersion and a systematic 
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component from observational uncertainties. These are 
added in quadrature: 



a{mk) 



rik 



( ily 



(17) 



where k = 1, ... ,17 labels the 17 bins, uo = 0.15 mag 
is the statistical calibrated uncertainty of an individual 
supernova, rik is the number of supernovae in bin k, and 
the redshift Zk is taken at the bin center, i.e. at Zk — 0.05 
for k = 1. The second, systematic term uses the error 
model discussed in [s^, with a linear rise with redshift. 
This serves as an approximation of observational uncer- 
tainties. SNAP instrumentation and observing strategy 
is specifically designed to limit systematic uncertainties 



to below nis 



0.02 mag out to the maximum redshift 



z = 1.7. 

The data is related to the theoretical parameters by 



m{z) = 5logi 



(l + z) / dz' {l-nD){l + z') 



/\3 



fine Jo 



3 /J"<^+^'' d{ln{l+z"))ll+wiz")]-] 



+ M 



While the use of only one constant parameter w to model 
the function w{z), i.e. assuming a constant equation of 
state, is physically unrevealing and often a poor approx- 
imation, more than two parameters makes it difficult to 
obtain meaningful constraints. We will write equation 
for m{z) using the fit Ijllf) containing two parameters, 
wq and Wa- 



m{z) = 51o| 
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(l + z) / dz' {l~nD){l + z') 
L 



/\3 



+ r!i5(l + z')'*'+"''+""'e~^'"°^ 



'1/2 



M . 



Note that the offset parameter Al = Af + 25 — 
51ogjQ(77o/100km/s/Mpc) (where M is the supernova 
absolute magnitude) just adds linearly, so the sensitivity 



derivative 



dM 



1 for all bins. Technically, other astro- 



physical terms enter into m{z) such as a change in su- 
pernova magnitude due to dimming by intervening dust, 
but we assume that we are dealing with fully calibrated 
data, with only the cosmological dependences remaining. 

In Fig. ^ we evaluate the sensitivity derivatives at the 
fiducial parameter values: flu = 0.72, wo — —1, Wa = 0, 
corresponding to the cosmological constant model for the 
dark energy. For the case of only supernova data we also 
add a gaussian prior of a^j^ = 0.03, corresponding to a 
reasonable level of knowledge on this parameter by the 
time SNAP data is available. Now we have all the ele- 
ments needed to generate the Fisher matrix. Next we in- 
vert it to obtain the covariance matrix and the parameter 
estimation uncertainties (the diagonal elements). If we 
want to plot the confidence contours in a two dimensional 



section of the parameter phase space, e.g. the wq — Wa 
plane, then we marginalize over the other parameters and 
reinvert the covariance matrix to get a reduced, 2x2 
Fisher matrix. This gives the major and minor axes of 
the elliptical contour (it is always an ellipse within the 
Fisher approximation) and the orientation, i.e. the direc- 
tion of the degeneracy between the parameters. This is 
basically an eigenvector where a particular combination 
of the parameters is best determined while an orthogonal 
combination is poorly constrained. 

A plotting program takes the eigenvector information 
and draws the ellipse, with the scale determined by what 
level of probability one wants to enclose within the con- 
tour. For 68% of the probability enclosed (called la 
joint probability), each axis of the figure is 1.52 (V2.30) 
times larger than the individual la probabilities (called 
la projected probability, or enclosing 39% joint probabil- 
ity) from the diagonal elements of the covariance matrix. 
For 95% confidence level, the scaling is 2.45 (\/5.99; for- 
mally the numbers 2.30 and 5.99 are the increments in 
the statistic from the fiducial model to a model on 
the 68%, resp. 95% contours). As found here, without 
a biasing systematic magnitude error, the best fit model 
will always be the fiducial, input model. 



Planck 

The cosmic microwave background data carries much 
information on the cosmological model parameters. Here 
we consider a simple subset of the measurements, in- 
volving only the distance to the photon last scattering 
surface, which is exquisitely determined by the location 
of the acoustic peaks in the CMB power spectrum. The 
Planck Surveyor mission will also precisely determine the 
combination VlMh^, so we employ the reduced distance 



dzf{z 



,-1/2 



where the upper boundary of the integral is taken at the 
last scattering surface (LSS), z^^g — 1089, and where 



(1 + zf + 



n, 



1 - n, 



3(l+«)o+tUa)g-3u)aTf7 



Planck has the ability to determine d up to 0.7%, i.e. 



CTj = 0.007 • d; this is translated into a Fisher matrix 



;iPlanck 



1 dd dd 
aj dp^ dpj 



(18) 



where p — {riu, "u^a}- Since we want to add this 
Fisher matrix to the Fisher matrix from SNAP, we will 
include also a fourth column and row of zeros (for the 
supernovae parameter M that does not enter the CMB 
data) . 
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The final Fisher matrix for SNAP and Planck com- 
bined then is 

^SNAP+Planck _ ^SNAP _|_ ^Planck ^^Q) 

Notice that it does not matter that the Fisher matri- 
ces are based on different data quantities from different 
experiments (m(z) and d). The information still adds. 
When including the Planck data, we can eliminate the 
step of adding a prior on flu, since the different cos- 
mological quantities measure sufficiently different combi- 
nations of fljy with the other parameters. This breaks 
degeneracies well enough to determine n d with superior 
precision (roughly equivalent to a prior of crjjj-, = 0.01 



Weak lensing 

The procedure for other data sets, such as from mea- 
surements of the gravitational distortion of images of 
background galaxies by foreground mass concentrations, 
known as WL, is implemented similarly. Here we used 
only estimates of the future precision of measurements of 
the linear part of the lensing shear power spectrum (see 
§3.3 of 34] for details). An even stronger data set em- 
ploying the full power spectrum, and possibly other lens- 
ing methods, should be available from the SNAP wide 
field survey and other experiments. 
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